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The claim that By 1 ^ — in Sec 4A of the above paper (hereafter as I) 
is incorrect. The calculation over there only establishes a particular solution 
Eqs. 1(33), 1(34) for a£} z to Eqs. 1(31), 1(32) with B y 1] = 0. It is necessary 
to add a solution to the homogeneous equations corresponding to Eqs. 1(31) 
and 1(32), i.e., 



dB 



(i),h 



y 



dB 



(i),h 



.7 



AW* = o (l) 

■AW' h =0 (2) 



correspondingly 



dx 

d 2 By 1) ' h d 2 By lhh (l)h 

-^ + ^-< ) ' /t = °- © 

This homogeneous solution must be included so that By 1 ^ — AirMy 1 ^ contains 
no discontinuity at the vortex. Using My 1 ^ = hence 

M v ] = , (4) 

the homogeneous solution required is 

B v ),H = -k^K^cose (5) 
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correspondingly By = By'* 1 . We see now that By 1 * 1 — AnMy 1 ^ becomes 
identically zero. To check the consistency with the Maxwell equation 1(8), 

we note that A^z' 11 can be obtained from Eqs ([1]) and Together with 
Eqs. 1(33) and 1(34), we obtained the correct first order vector potentials as 

I« = Kp-K (r) (6) 

Z7T 

and 

A[ 1] = (7) 

Eq © gives J { z 1] = 0. Eq ©, together with 1(1) and B y 0) = ^-K Q {r) gives 
4 1} = 0. 

The total field up to first order is therefore 

4 0) + B v ] = % \ K ^ r ) - S*i(r)cos0] (8) 
k ^.[K {\v + kz\)] (9) 

Therefore, to first order in k, the field pattern is the same as a cen- 
trosymmetric superconductor except that it is shifted from the vortex center 
towards — z by k in dimensionless units (k\ x in ordinary units). This shift 
has been found numerically by Oka et al pQ in their numerical solution of 
Ginzburg-Landau equations, instead of the London equations here. 

In contrast, the total gauge invariant superfluid velocities along x and 
z are, up to first order in k, proportional to —Ki(r)cos9 + kKo(r) and 
Ki(r)sin9, respectively. The singularity is dominated by K\(r) terms and 
is still at the vortex center. 

The correct Ax/z in Eq © and has to be used to obtain the next 

(2) 

higher correction B y , which can be shown to be finite instead of zero as 
claimed in I. The expression is not particularly informative and we shall not 
obtain it here. 

The results of the other sections in I are unaffected. 
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Abstract 

We study the magnetic properties of a superconductor in a crystal without z — > — z symmetry, 
in particular how the lack of this symmetry exhibits itself. We show that, though the penetration 
depth itself shows no such effect, for suitable orientation of magnetic field, there is a magnetic field 
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discontinuity at the interface which shows this absence of symmetry. The magnetic field profile of a 
vortex in the x — y plane is shown to be identical to that of an ordinary anisotropic superconductor 
to second order in a small parameter k. For a vortex along z, there is an induced magnetization 
along the radial direction 
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I. INTRODUCTION 



Lately, there is much attention on the theory of superconductivity in systems without 
inversion symmetry in the normal state, e.g.,-^^^ This surge of interest is due in no small 
part to the discovery of new superconductors in crystal where this symmetry is absent, for 
example, the heavy fermion superconductor CePtsSi 8 with T c ~ 0.75K. The normal state 
of CePtsSi (ignoring the possible antiferromagnetic ordering at T/v ~ 2.2K) has tetragonal 
PAmm symmetry. Due to the displacement of Si along c direction, the reflection symmetry 
z — > — z has already been already lost in the normal state before the superconducting 
transition. 

Indeed on very general grounds, the properties of a superconductor in crystals without 
inversion symmetry is expected to be very different from those where such a symmetry is 
respected. In the latter case, which is applicable to most known superconductors, inversion 
symmetry and Kramers degeneracy allow the classifications of superconducting states into 
singlet, even parity on the one hand and triplet, odd parity on the other 9 . The physical 
properties of such superconductors can then be drawn from the generalization of theories of 
conventional superconductors or superfluid 3 He. This is no longer the case if inversion sym- 
metry is already broken in the normal stated Some peculiar aspects of these superconductors 
have already been discussed theoretically in the literature.— 

In particular, in ref^i, it is shown that for systems where z — > —z is broken, a superfluid 
flow in the plane, say along x, is expected to produce a spin polarization along y when a 
Rashbaii type spin-orbit coupling is present. Conversely, a Zeeman magnetic field along y is 
expected to generate a superfluid flow or phase gradient along x. Though the calculation in^ 
was specifically for two dimensions, the general argument there is obviously generalizable to 
a three dimensional superconductor, (see also Section |V] below) Here we elaborate further 
on the consequence these effects discussed in Ref^. We shall (Section IIIip first study the 
screening of an external magnetic field by the superconductor, i.e., the Meissner effect. More 
specifically, consider a magnetic field along y with the superconductor occupying either z > 
or z < 0. These two situations are not equivalent if the crystal lacks the z — > —z symmetry. 
We shall however see that, (at least to the surprise to the present author), the penetration 
depth itself shows no direct effect of the absence of this z — > —z symmetry. The latter only 
manifests itself as discontinuities in the magnetic field near the crystal surface with opposite 



2 



signs for the two mentioned geometries. 

Next (Section HVj) we discuss the structure of a vortex in London theory. We shall show 
that, at least to second order in a small parameter k to be defined below, a vortex for 
field along y has a magnetic field distribution again symmetric with respect to z —>■ —z. 
Hence the constant magnetic field contour is elliptical similar to that of an ordinary crystal 
with different effective masses along x and z. For magnetic field along z, there is a spin 
magnetization pointing along the radial direction, the sign of which reflects the absence of 
z — > —z symmetry. We finally estimate the order of magnitude of these broken symmetry 
effects (Section IVj). 



II. CONSTITUTIVE EQUATIONS 

First we recall the constitutive equations^ relating the (number) current J, (local) mag- 
netization M, gauge invariant phase gradient q = KV(p + —A (electron charge = — e) and 
the magnetic field B in a superconductor with broken z — ► —z symmetry in its normal state. 
For currents and magnetic field in the x — y plane, the relations are expected to have the 
general form 



J x = n s — kB v (1) 

2m x y V ; 

J = n s -^- + kB x (2) 
y s 2m y K ' 

Hi 

M x = ~-q y + XxB x (3) 
M y = \q x + XyB y (4) 

Here n s is the superfluid density and m x etc are the effective masses, and \x etc. the 
spin susceptibilities (m x = m y and Xx = Xy)- The terms proportional to k in eq ([H) and 
d2J) represent the supercurrent induced by Zeeman field^ and the corresponding ones in eq 
([3]) and (|3J) represent the magnetization induced by superflow^. Note the difference in sign 
for the terms proportional to k between eq ([T]) and (j2j) and between eq ([3]) and (jlj). These 
terms are specific to the absence of z — > —z symmetry. The purpose of the present paper 
is to study effects due to these terms. As already mentioned in ref-, these equations are 
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expected from the general form of the free energy 



im x am y sm 2 

\x,Bl + \ X ,Bl + \ X ,Bl) + ^ (5) 
~2 (<lxBy - q y B x ) 

appropriate to the present symmetry . The term proportional to k here is symmetry allowed 
in the present case (see also^ 1 ^). 

Currents and magnetization along the z axis are given by the usual relations 



M z = Xz B z (7) 

The equation governing the magneto-statics is given by 

V x B = 4ttV x M + -y(-e)J (8) 

where J and M are given by eqs (CD)- (jlj), (E]) and (JTj). Eq (jSJ) also follows from the free 
energy (jSJ) by variation with respect to A when one takes into account the basic relation 
B = V x A. 

It is worth mentioning that the screening of a magnetic field in a superconductor without 
inversion symmetry has also been considered in ref^. However, they assume the presence 
of a term with J || B. This term is absent for our case of P4mm symmetry because the 
presence of vertical reflection planes such and y — z and the fact that J and B 

transform differently under reflections. 

III. MEISSNER SCREENING 

(A) Let us now consider screening of a magnetic field in the basal plane, taken to be B = B y y 
without loss of generality, with the sample occupying either z > or z < 0. These two (not 
a priori equivalent) geometries are particular interesting since the broken z — > —z symmetry 
may manifest. With translational invariance along x and y, one can verify that the z 
component of eq (jHI) is trivially satisfied. The y component is satisfied by B x = 0, M x = 0, 
J y = 0. The x component reduces to 

^ = ^ + *n Jm (9) 

oz oz c 
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In our case we can take the gauge = and A = A x (z)x. With eq p]) and (TjJ and taking 
into account B y = it can be seen that the terms proportional to k cancel out in eq ([9]). 
Further taking the z derivative gives 

d 2 B y _ Aixn s e 2 
dz 2 m T c 



4 ^)^T = — V 5 * • ( 10 ) 



Thus the penetration depth A^ is given by 

1/^ = -^/(l - 47TXV) (11) 

Here the subscript of A x denotes that the current is along x. Typically ^ 1 and thus \ x 
reduces to the usual expression l/X 2 , = ■ 

Since the term linear in k drops out, the penetration depth shows no direct effect of the 
lack of z — >• —z symmetry. It is the same for samples occupying z > or z < 0. 

This, however, does not mean that there is no broken symmetry effects at all. Since 
By — AirMy has to be continuous across the vaccum-sample interface, we then have, for 
sample occupying z > 0, B cxt = B m — 4irM y (0 + ). Here B m = B(z = 0+) is the value of 
magnetic field just inside the sample. M y is given by eq (EQ), hence (recall that = 0) 

My(0 + ) = ^A x (0 + ) + Xy B m (12) 

Inside the sample, B y (z) = B in e~ z ^ Xx with A x already determined in eq (TTTT) above. Thus 
A x (z) = -\ x B in e- z / x *. We finally have 

B in = B«t/(1 - 47TX, + R) (13) 

where R = ^f-\ x is a dimensionless parameter. 

Similar calculation for the case where the sample occupies z < shows 

B in = B cxt /(l-4n Xy -R) (14) 

Hence the k term results in a discontinuity of the magnetic field with a contribution of 
opposite signs in the two geometries. This is a manifestation of the broken symmetry. 
Of course the magnetic field is discontinous only under our (London) approximation: the 
variation would spread out probably over a length scale of order of the coherence length in 
a more microscopic treatment. 



5 



(B) For ease of latter reference we also consider the screening of magnetic field B = B y y 
with the sample occupying either x > or x < 0. Since there cannot be any current along 
x, J x = and so from eq (JTJ) 

n q e ( . hcd(b\ ^ 

A x + -Z)= *B y (15) 



m x c \ 2e dx y 

Note it follows that we cannot set both A x and |^ zero, in contrast to the more usual case 
where the k term is absent. In our case, it is convenient to make use of the translational 
invariance along z and choose A x and 4> both dependent only on x, so that, e.g., B y = — 
Eq (jSJ) gives 

f, = 4 ™,_^ (16) 

ox ox c 

Eq (TJJ, together with eq (ITB1) . gives 

My = (Xy + ^)B y (17) 

Therefore eq f|T6|) reduces to 



(1-4^-^ = -^^ (18) 
Taking the z derivative shows that the penetration depth A 2 is given by 



1/^ = ^/(1-4^-^) (19) 



Thus the k term only gives a correction to the penetration depth proportional to k 2 . Hence 
again there is no asymmetry between the geometries where the samples occupy x > or 
x < 0. Similar argument as in the last subsection shows that there is a discontinuity in 
magnetic field near the sample surface oc k 2 . 

(C) We finally consider a field along z. Without loss in generality, we take the sample 
to occupy y > 0. Translational invariance along x and z are respected and all quantities 
depend only on y. It can be shown easily that the magnetic field obeys the usual screening 
equations and thus B z {y) = B z (0)e~ y ^ Xx where 1/X 2 = — ^Xz) is the same as 

that in section (A) [except that the small correction due to spin susceptibility is here now 
(1 — 4:7iXz) instead of (1 — ^Xy) m ec l (flTj)] . The peculiar feature here, however, is that 
from eq (0} that M y ^ since q x ^ 0. One easily finds M y {y) = ^B z (0)e~ y/Xx . Thus there 
is a magnetization towards (if k > 0) or away from (if k < 0) the inside of the sample if 

-Bzext > 0. 



IV. FIELD DISTRIBUTION OF A SINGLE VORTEX 



(A) Now we study the magnetic field profile around a vortex for magnetic field in the 
basal plane, chosen to be again along y. For simplicity, we shall ignore the small spin 
susceptibilities Xy m the present section. The basic equations are again eq ([HD and (fTET) . 
with My, J x , J z given by eq (|4j) , (jTD , ([6]) . In constrast to section UTTlA the vector potential 
must depend both on x and z, and thus 

B ^-^x &A Z f 2()1 

v dz dx 

Due to the presence of these two contributions to B y , the terms linear in k does not drop 
out in eq (Q (or equivalently eq fl26l) below). Therefore it is not a priori obvious that the 
vortex field distribution will obey z — > —z symmetry. We shall however show below that 
this symmetry is respected at least to order k 2 . 

We begin by performing a rescaling of coordinates by the penetration depths, thus we 
write 

z = \ x z (21) 
x = X z x (22) 

It is convenient also to rescale the components of A: 



A x — X X A X (23) 
A z = X Z A Z (24) 



so that 



dA x dA z , . 

B -lF-aT (25) 

In these variables, eq ([9]) and eq (fT6j) become 

where $o = ^xiT ^ s a sca l e< i ^ nx quanta (magnetic field). Note that as usual, B and J 
must vanish at large distances and the total flux is therefore given by the flux quanta 
In our scaled variables this condition becomes J B dxdz = $ . 
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It seems difficult to solve eq (|26|) and fl27j) for general k. We shall thus make use of the 
smallness of k to solve these equations order by order in this parameter. (Strictly speaking 
the penetration depth A 2 already contains a k 2 correction, but the present rescaling using 
this corrected A 2 simplifies the calculations below substantially). 

In the lowest (zeroth) order, we can drop all terms on the right hand sides of eq (I26p and 
(|27|) . The resulting equations are the standard equations for the vortex. For a vortex with 
positive flux along y, we need to choose = —9 and the solutions are 



where 



A 



(0) 



a<°> = Afe 



— ~-Kx(r) 
2n\r u J 



(28) 
(29) 

(30) 



and we have introduced the cylinderical coordinates for the scaled variables: r = [z 2 + x 2 ) 1 ^ 2 , 
cos# = z/r, sin9 = x/r, 9 = — sm9z + cos9x. K , K\ are the modified Bessel functions. 
To first order, we have, from eq (I26I) and (T5j 



(31) 
(32) 



dB v 1] _ _ z.^_ ( 2(0) , £o^ | _ ~ n an 
dz K ~ K d~z \ A * + 2*8x1 y 



dB 



(i) 



dx 



- k l. f i(o) + £^ 

A - ~ dx \ Ax + 2ndx 



We claim that By = 0. Assuming this, using eq (EE}, (EHD and therefore + fa 
^Ki(r)cos#, eq (J3TJ> and ([32]) become 



(0) , $o 



()'i 



27T 



K (r) + |r(^i(r)cos^) 
az 



2tt 



9 

— (i^i(r)cos6 



(33) 
(34) 



Indeed, substituting these equations into (1201) . we get 



9 2 



C) 2 



<9z 2 9a; 2 



(#i(r)cos0) 



(35) 



However, by properties of the Bessel function, ( ^ + ) (i^i(r)cos^) = K\(r)cos9. Fur- 



ther using = —Ki shows that indeed B 



(i) 



' \dz 2 
0. 
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Now we proceed to the second order. Eq ([26]) and ([27)1 read 



dB 



(2) 



dz 



4(2) _ 

92; 



, 1(2) . . -li(i) 



ft 



1(0) , fo^ 
A « + 2vr<95 



(36) 
(37) 



We again claim that By 
and flS}, 



(2) 



0. If so, we get, using Af ] + gff = |^i(r)sm#, and eq ([33 



<J>0 



-ft 2 



2tt 



if = ft 2 ^ 
2tt 



dK (r) ^ s6 
dx dxdz 1 



From eq ([20]) for 5^ ; , we get 



y Z7T 



V 2 K (r) + ^rV 2 (iT 1 (r)cos6 



lfi(r)sin0 



.9 

— (Ki(r)sin6 



(38) 
(39) 



(40) 



Using again the properties of the modified Bessel functions, this reduces to 

B^ : 



-2*0 
-ft 

2tt 



K (r) + — (#i(r)cos0) + ^ (Kx(r)sin6 



(41) 



The last two term combines to ^ 4- (rKx(r)) = —K by recursion relation of modified Bessel 



function. Hence B 



(2) 



as claimed. 



Hence to second order in ft, B y has the same form as an ordinary anisotropic material: 
By = 2nX >0 x K ([(z/\ x ) 2 + (x/Xz) 2 ] 1 ^ 2 J . The constant magnetic field contours are ellipses 
with center at the point where the order parameter has a singularity. 

We would however like to add two cautionary remarks. Firstly, it is not true that other 
physical quantities such as J or M also have elliptic distributions around the vortex. For 
example, using eq (j3J) and the solution to A^°\ we find My 1 ^ = —k-^Ki(r)cos6. Therefore 

—z. Secondly, it is not true that the corrections to magnetic 



M y (x,z) is odd under z 



(3) 

field vanish to higher orders. Proceeding to the third order, one can verify that B y ' = 
is inconsistent with A( 2 ) of eq (1311 and (13"8"|) . We, however would not proceed to calculate 
these very small corrections. 

(B) We now consider a vortex with field along z. Using translational invariance along z, one 
finds that there is no magnetic field induced along the x — y plane, and B z and A are the 
same as those of an ordinary superconductor with in-plane penetration depth X x given in Sec 



IIHI before. However, due to the presence of q, there is an in-plane magnetization induced by 
the presence of the k term (see eq ([3]) and (j3])). We find that this magnetization is along the 



magnitude of this magnetization is therefore k/Air times the local magetic field B z along 
the z direction. It points radially outwards if n > 0. For field along —z, this magnetization 
changes sign and points radially inwards if k > 0. 

V. ORDER OF MAGNITUDE 

Finally we estimate k, assuming the clean limit. The calculations in Ref^ can be easily 
generalized to the present 3D case once the Fermi surface and the spin-orbit splittings are 
given. We however would not do this calculation here but just satisfy ourselves with some 
estimates. It should be noted that, for the crystal symmetry PAmm in question, the allowed 
spin-orbit interaction, besides one in the Rashba form —az ■ (p x a) = —a(p x a y — p y u x ), 
(here p is the momentum, a the Pauli matrix, and a a coefficient) can also have terms of the 
form — /3p x p y {p x a x — p y a y ) and —JPxPyPzipl —pfy&z- Here a, j3, 7 can be functions of p but 
must be invariant under the symmetries of the crystal, (c.f., 13 ). Both the first two terms can 
generate the terms proportional to k in eq (Jl])-(j4j). However, we expect that (c.f., eq (12) of 
re&) that k (at T — > 0) has an order of magnitude given by ~ Hb(Pf+ ~Pf-)/^ 3 where /ib 
is the Bohr magneton, pp± are the typical Fermi momenta for the spin-orbit splitted Fermi 
surfaces. Using the expression for and n s ~ p 3 F , we find that k is of order 



where 5 is the typical splitting in energy by the spin-orbit interaction, /z the chemical po- 
tential, and tb the Bohr radius. For 5/fi ~ 0.1, this ratio is then of order ~ 10 -3 assuming 
typical electron densities. 

VI. CONCLUSION 

In conclusion, we have studied some magnetic properties for a superconductor with no 
inversion symmetry in its normal state. In particular we investigated how the broken symme- 
try and magneto-electric effects discussed in ref^ 1 ^ exhibit themselves in Meissner screening 



radial direction: M = M r f (thus its curl vanishes) where M r = k 



^(^).The 




(42) 
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and vortices. An unusual magnetization spatial pattern are found in some geometries. This 
magnetization can in principle be detected by Knight shift measurements. 
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